Source of Acquisition 
NASA Marshall Space Flight Center 


Chaotic motion of relativistic electrons driven by 
whistler waves 

G V Khazanov 1 , A A Tel’nikhin 2 , and Tatiana K. Kronberg 2 

1 NASA Marshall Space Flight Center, National Space Science & Technology Center, 
320 Sparkman Drive, Huntsville, Alabama 35805, USA 

2 Altai State University, Dimitrov, 66, Barnaul 656049, Russia 
2 Altai State University, Dimitrov, 66, Barnaul 656049, Russia 

E-mail: tatiana . kronbergSlist . ru 

PACS numbers: 52.65, Cc 52.20,Dq 05.45,Pq, 94.20,Rr 



Chaotic motion of relativistic electrons driven by whistler waves 


2 


Abstract. Canonical equations governing an electron motion in electromagnetic field 
of the whistler mode waves propagating along the direction of an ambient magnetic 
field are derived. The physical processes on which the equations of motion are based 
arc identified. It is shown that relativistic electrons interacting with these fields 
demonstrate chaotic motion, which is accompanied by the particle stochastic heating 
and significant pitch angle diffusion. Evolution of distribution functions is described 
by the Fokker-Planck-Kolmogorov equations. It is shown that the whistler mode waves 
could provide a viable mechanism for stochastic energization of electrons with energies 
up to 50 MeV in the Jovian magnetosphere. 


1. Introduction 

A description of the stochastic dynamics of charged particle in the field of wave packet 
is one of the fundamental problem in the theory of plasma physics, that has attracted 
current interest for many years [1]. Thus the stochastic dynamics of relativistic electrons 
in the time-like wave packet has been discussed by Chernikov et al [2]. A theory of 
dynamic chaos has also been employed to deal with high-frequency heating of electrons 
in the Langmuir wave packet [3]. Possible astrophysical applications of this results were 
discussed in [4], Nagornykh et al have been announced some results of theoretical study 
concerning stochastic particle motion in the magnetized plasmas [5]. On the other hand, 
Zaslavsky et al [6] have been reported the results describing electron stochastic motion 
in the upper-hybrid wave propagating across an ambient magnetic field. Concepts 
and methods which were developed in these works have impacted our approach to the 
problem. 

The goal of present paper is a systematic study of relativistic particle motion in a 
whistler wave packet. Earlier the basic concepts of energy diffusion of relativistic 
electrons resulting from resonant interaction with whistlers in the magnetosphere have 
been discussed, for example, by Walker [7], and Summers et al [8]. In our work we discuss 
stochastic motion of high-energy electrons due to a spectrum of whistler mode waves. 
The results of study of the problem are of interest from the point of view interpreting 
the experimental results. 

The paper is organized as follows. The canonical equation of motion expressed 
in terms of the action- angle variables are derived in Sec. 2. The nonlinear resonance 
wave-particle interaction is discussed qualitatively in Sec.3. In Sec.4 it is shown that the 
stochastic motion of relativistic electrons can be described by the closed set of nonlinear 
difference equations. The solutions of these equations are obtained both analytically 
and numerically. In Sec. 5 it is proved that particle dynamics is realized on stochastic 
attractor with fractal structure. The Fokker-Planck-Kolmogorov equation is derived in 
Sec. 6, and the effects associated with particle stochastic heating such that the pitch angle 
scattering and radial drift of plasma particles are studied. Application of our results 
to high-energy electrons observed in the Jovian election radiation belts is described in 
Sec. 7. In Sec. 8 we give the conclusions of our studies. 
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2. Basic equations 


Let us consider a relativistic particle of charge \e\ and mass m in the wave packet of 
extraordinary electromagnetic waves propagating along an external uniform magnetic 
field of strength B. The Hamiltonian corresponding to the problem is 


H( r, p; t) = \jm 2 + (p + A) 2 , 

(1) 

and the canonical equations of motion are 
p = [p,i7], r = [r ,H], 

(2) 


where p is the canonical momentum, r is the position vector, A the vector potential 
divided into two parts, A = A w + A ext , superscripts w and ext denote both the wave 
and external fields, and [ , ] stand for the Poisson brackets. 

We have employed here and throughout this paper, the frame of reference in which the 
speed of light c = 1 and charge \e\ = 1. 

We denote by R the set of all real numbers. Then p € R 3 , r € R 3 , and the smooth 
manifold M — R & will be a canonical space of this dynamic system, and Rf = R 3 x R? 
is a direct product space. 

In order to write down equations of the particle motion one must specify a 
coordinate system. We have chosen a Cartesian spatial coordinates system whose 2 axis 
is directed along the external magnetic field, the plane perpendicular to this direction 
is spanned by the orthogonal coordinates x and y. 

Making use of the connecting relations, 

B = [V, A], E = -dA/dt, (3) 

we have in the coordinate representation 


r =(x,y,z), B ext = (0, 0, B), 

A w = [Y, A k sin<y9, ]TH fc cos<£,0 j , 

\ k k ) 

c p — zk — talk, 

A ext = (-By, Bx, 0)/2. 


( 4 ) 

( 5 ) 

( 6 ) 


Here the expression for A ex£ is written in the axial gauge, A * is the amplitude of mode 
in the wave packet, k is the wave number, and w* is the dispersion equation. 

The dispersion relation for electron branch of the whistler mode waves in the cold 
magnetoplasma is written as 

k 2 /u? = 1 + Wp/[w(w B - w)], (7) 


where and u p are the gyrofrequency and electron plasma frequency, respectively. 
Equation (7) in a weak magnetic field when the condition (w^w/cu 2 ) <C 1 is valid, reduces 
to 


V = 


w(wg — w) 




v ph < !• 


( 8 ) 
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In a strong magnetic field, provided that uj/ujb — oj p /o>b, Wp/(u>u; B ) <C 1, equation (7) 
reduces to 


Vph 1 



2(jJU) B 


( 9 ) 


We now take into account axial symmetry of the non-perturbative problem, and 
introduce the new variables, an action(/), and an angle (0), by a canonical 
transformation (x,p x \y,p y ) — > (9,1) : 


x=rcos9, p x = — (mrwj3/2) sin 0, 

y — r sin 9, p y — (mrco B / 2) cos 9; 

r — \j2muiBl /rnojB,oj B = B/m, 
where r is the gyroradius. 

The Hamiltonian (1) in this representation becomes 

H(z,p z ; 9,1; t) = H 0 (p, /) + \J‘lmw B IHf ) x 


(10) 

( 11 ) 


' cos(zk + 9 — twfc), (12) 

k 

H 0 (p,I) = i/m 2 +p 2 + 2 mu B I, Pz^P- (13) 

Here we have assumed that the ratio p — A w jm <C 1 is the small parameter of the 
problem, and retain in (12) only the leading terms. 

Associated with (12) the equations of motion are 

p = [p,H] = \j2rmz B IH^ -^kAk sin ip, (14) 

k 

i = [I, H] = yj2 muj B I H q 1 - A k sin ip, (15) 

k 

z = [z, H] = pH^ 1 , e = [0,H] = u> B mH. (16) 
We omit in (16) the terms of the order of // 2 and introduce the definition for the phase 
ip d = zk + 9 — u) k t. (17) 


It is obvious that a Hamiltonian flow (15, 16) preserves the measure dp, = dp z dldzd9 
on the smooth manifold M = B? x S, where S is the circle and 0(mod27r) 6 S. Two 
pairs of the canonical variables p, z and 1, 9 are a symplectic coordinates on M with 
that algebra: 

[fipi = 6/] =: [ e,j,e z ] — [ez>e#] — 0, 

\e p , e z J [6/760] 1 5 


where [,] denote a skew-product, the set ( ej ) is the coordinate basic, which can be given 
by a differential operators d/dx P (Schutz, 1982). 

This algebra will be used in the following. 
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Now the structure of wave packet, A w (t,z) — Jfk exp[i(zk + 0 — tu>k )\ , is to be 
specialized. The most frequently used representations of wave packet are the so-called 
time- and space-like representations [6], [!)]. Such wave packet may be exited in a plasma 
due to the intrinsic instabilities. In this particular physical situation, we make following 
simplifying assumptions regarding the structure of the packet, 

k = k 0 + nAk, u = w 0 + nAw, Ak/k 0 <ttl, Aw/w 0 < 1, n € Z, (18) 

where ko,u>o are the characteristic wave number and frequency, Aw(A/c) is the group 
dispersion in wave spectrum, so that 

Aw = v gr Ak, Vg r —du/dk, (19) 

v gr is the group velocity, 

Ak — 27 r/L, A to — 2ir/T, (20) 

L, T are the length- and the time- scales of the problem, Z denotes the set of all integers. 
Then we suppose that the characteristic spectral amplitude A 0 is a slowly varying 
function on t and z such that 

A 0 /u)A 0 ^ 1/wT (< 1), X7A 0 /kA 0 ~ 1/kL (< 1), (21) 

and write down the wave packet in the form 

A w (t, z) = Aq exp(i^o) Y Ak ex\>[i{7iAkz -I- nAvot)\, (22) 

ti£Z 

tp 0 — k 0 z + — u)o)t. (23) 

Define the parameter, rj namely, the ratio of the particle velocity along an ambient 
magnetic field to the group velocity, 

r\ = v z /v gr . (24) 

In the limit, Ak -> 0, — > 0, expression (22) can be transformed into 

A w (t, z) = 2 i 0 exp(iV’o) ~ n ) ( 25 ) 

TldiZ 

There is the time-like representation (TLR) of a wide wave packet. The Poisson sum 
formula 

Y. exp(inAcot) = EV-") (26) 

nGZ nEZ 

has employed in (25), where 5(-) is the Dirac delta-function. 

In another limit Aw — > 0, 7] -> oo, we can easily show that the wave field takes the 
form of the space-like (SL) wave packet 

A w (t, z) — yl 0 exp(#o) n )• ( 27 ) 

nGiZ 

This packet represents a periodic sequence of impulses with characteristic spatial period 
L — 2ir/Ak. 
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Note that both representation are often used. Thus it is established [6], [9], that the 
TLR is available for the problem if the condition r q 2 - Cl holds. The TLR of the electric 
field of the electrostatic waves was used by Chernikov et al [2] to derive the relativistic 
generalization of the standard map. On the other hand, the SLR has been utilized in 
[3] to describe the stochastic motion of relativistic particle in the electrostatic field of 
Langmuir waves, whose group velocity is small as known. 

3. Qualitative analysis of the particle motion 

We now proceed to the physical aspect of the problem. First we consider the dynamics 
of an autonomous system with the Hamiltonian given by (13). In this commutative 
integrable situation (1 = 0, p = 0), the phase curves lie on constant energy surfaces of 
Hq = E — const: 

{(z,p;9,I) : \Jm 2 + p 2 + 2 mco B I = E}. 

All solutions are regular, since trajectories in phase space cannot intersect and each 
phase curve correspond in configuration space to particle gyration in an external 
magnetic field with a frequency 

To understand the physical picture of the stochastic motion, let us discuss the simple 
case of the single-wave Hamiltonian, 

H = H 0 (p, I) 4- \j2mu B IHf l ■ A k cos (zk + 9 - u> k t), (28) 

where A k is the amplitude of certain mode in the wave packet, k is an arbitrary 


wavenumber, u k is the frequency of this mode. 
Thus the associated equations of motion are 


p = k^muBlHf 1 ■ A*;Sinr />, 

(29) 

/ = \j2mu>BlHC l ‘ Ah sin-0, 

(30) 

z = v z = pH 0 1 , 6 = 

(31) 

'tph = w(fc, E) = kv z + losttiE -1 — oj k , 

(32) 


<jj(k, E ) is the nonlinear phase frequency. 

It is easy to see that the system possesses the first integral, inv. = (p — hi) = 0. We 
put const. = 0 below and represent the invariant of motion as 


ap — 2 oobI = 0, a = v P h- (33) 

to 

It is well-known [1] that the stochastic motion will be mainly determined by the 
nonlinear resonance wave-particle interaction, 

u)(k, E r ) = ( kv z + cobtuE ~ l — Lo k ) r = 0, (34) 

as well as the group velocity dispersion. 

In the wave spectrum case the equation given above determines the family of resonance 
states in energy space. The phase space of resonance states has structure of the product 
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bundle, BxF , where B is the base, i.e. the open subset U C i? 2 , I r ,p r £ U, and F is 
the typical fibre, F = T 2 , with the coordinates ((u> — kv z )t, #(modd 2tt)) G T 2 . 

Now varying of (34) with respect to u, k, E r we obtain the interval between an adjacent 
fibres, 

A E — (v z A k — A u)/(h(dv z /dE) r — u B m£" 2 ). (35) 

Here A k and A w are the intervals between an adjacent modes in the wave spectrum, 
and A E is the interval between the two adjacent resonance states in the particle energy 
space. 

We specify the problem. First we introduce the following notation 

e = E/m, £ z =p/m , e t = ^2mu B I /m, (36) 

and denote by sup{e} = the upper bound of the set {e}. 

On condition that (8) is valid, and the parameter a is small, we stipulate 

el max {a 2 , 1}; (37) 

unlike this case, we assume that the condition 

1< el < a 2 , (38) 

suits to the problem, when a 1, and (8) holds. 

Now with respect to the wave spectrum. The term w d = A kv z — Aw occuring in (35), 
describes the group velocity dispersion, and it can be represented in the form 

u d — A kv z ~ L~ l v z , (39) 

for v z /v gr 1; and, 

u d = Aw ~ T _1 , (40) 

for v z /v gr <C 1. 

Note that these representations are tantamount to the SLR and TLR of wave packet, 
respectively. 

We are interested in how relativistic and stochastic effects modify the property of particle 
acceleration in the case corresponding to particle motion in a SL wave packet. 

In the ultra-relativistic limit, both conditions, (37) and (39), are valid, and v z ~ 1, 
consequently, from (35) the following expression results 

Ae~e 2 (w bL)~ 1 , (41) 

in which the relation A k ~ 1 /L has used. 

Resonance wave-particle interaction leads to a widening of the own resonance energy 
level. Denote by 8e = e — e r the width of the resonance level. 

Integrating equations of motion (30) and (31) with the help of invariant of motion (33), 
we estimate the level width 

Se ~ ka 1 ^ 2 ef 1 ^ 2 (A k /m)At, 

A t=L/v z , v z ~l. 


(42) 
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Now the criterion, defining the allowed domain of random motion in energetic space can 
be written as 5e > As [10], [9] or, in the explicit form, 

£ <e b - (a 3 * ^(kL) 2 (A k /m)/2y /5 , (43) 

which follows from the two above results, (41) and (42). Formula (43) predicts that the 
upper value of energetic spectrum depends on the magnitude of wave field as A 2 / 5 . 

For describing stochastic behavior of relativistic electron in a TL wave packet, 
we use the conditions (8), (40) and (35), which correspond to the particular physical 
situation. Now substituting (40) into (35) we find the interval between the two adjacent 
resonance states, 

A £ ~ e 2 (argil) -1 . (44) 

Then with the help of equations of motion (31) and (32) we evaluate the level width 

5s ~ Lo B T(A k /m)/2s r . (45) 

Now the condition of resonance overlap, 5s > As yields the upper bound of the energetic 
spectrum, 

s b = (a 2 u 2 T 2 A k /7nj /3 . (46) 

Our analysis infers the possibility of self-generated, chaotic behavior of the system in 
some region of energetic space. It is easy to see that the structure of resonance energy 
space mentioned above corresponds to quantized energy space, when each resonance 
state is assigned the own energy level. The interval between adjacent energy levels 
caused by the group velocity dispersion depends on the particle energy as Ae oc s 2 . 
On the other hand, the width of each energy level is governed by the resonance wave- 
particle interaction, which becomes weak at very high energies of the particle. Thus the 
overlap criterion determines the allowed stochastic heating domain. We will carry on 
our advanced studies of stochastic motion in the next sections. 

4. Relativistic particle dynamics in a coherent wave packet 

We specify first the wave spectrum of a packet. Let us assume that the wave packet 
is given by (27). Thus, we will consider a relativistic electron motion in the space-like 
packet (SLP) of the whistler mode waves. In this approach dropping the subscript ”0” 
we write down the equations of motion (15) and (16) in the form: 


p = kA\j2mio B IH 0 'sinip^SiC-n), (47) 

nSZ 

I = \J 2 rnuj 11 IHf 1 A sin ?/; E *K - n). (48) 

nez 

z =pH ( f 1 , 0 = ojsmHff 1 -, (49) 

H 0 (jp, I) = \jm 2 +p' 2 + 2moj B I, (50) 

if = u(p, I) — kpH^ 1 + ojBmHpf 1 — oj. (51) 
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Here A, to and k are the magnitude, frequency, and wave number of the fundamental 
(characteristic) mode, ( — (z/L), L is the characteristic spacescale, S n = <5(£ — n), d(-) 
is the Dirac delta function, and Z denotes the set of all integers. 

The evident symmetries of rotation, if — > ?/ ; exp(iV’o)) Cq is an arbitrary constant, 
and a translation, {p, /}(C) -> {p, /}(C+1), allow us to represent dynamics of the system 
as certain iterative process by identification of the planes n£ and (n + 1)((. The group 
symmetry is realized first as the invariant of motion, 

p — kl = inv. (52) 


By means of (36), we write down expression (52) as 
- oT l e] = 0, a= (2u B /to)v ph , 


(53) 


where the constant of integration has chosen equals zero. In view of (53) the number of 
dimensions reduces to two. Let the variables z and p will be the represented pair. 

Now the explicit form of iteration system is to be found. Denote by TM the tangent 
fibering. There exits the map v : M — > TM, such that the field of vectors has the form 

d_ 
dz' 


v 




+ pHo 1 ir.- 


(54) 


Then each element of v is assigned the 1-form, so that the field of 1-forms may be written 
as 


v = -k^^Asin V 5 n dz + pH 0 'dp. (55) 

no 

The notations d/dx l and dx 1 are employed here for the orthornormal coordinate basis 
of vectors with algebra mentioned above, and for the dual basis of 1-forms, such that 


9 _ ,j, _ [ 1, i = i 
dx i 1 \ 0, i 7^ i' 


where (a;*) = (p, z). 

It is well-known that section of the field of 1-forms, v(v) = 0, on a submanifold in M is 
equivalent to solution of equations (48, 49) ( Schutz, 1982). 

Then direct integrating the equation (—ky/2mujBl H q 1 A sin ip 2 8 n dz + pHf l dp) (v) = 0 
along with (53) yields 

(p n+ i/m) 3/2 = ( Pn/mf 1 2 + (3/2)a 1/2 N(A/m) sin -0^, N = [ kL ], (56) 


where N is the characteristic number of modes in the wave packet, [ ] denotes the integer 
parts. 

Taking into account that the timescale of the problem is T — Ljv z , from (51) follows 
the equation bVi+i = b’n + io(p n+ i, I n \ \ )T, which describes the phase shift acquired by 
particle. This equation may be written in the explicit form as 

</Wl = <!>n+N (l + ^ (l - ^\/ 1+ ( ? ^ 1 ) 2 + “ ? ^ i )) < m ° d 2?f )'< 57 ) 
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It is readily seen, if the inequality e 2 1 is valid, the set of equations (56) and (57) 
can be transformed into the map: 

u n+ i =u n + Q sin tp n , 


•*/Wi = 077 + (37r 5/3 /2Q)|u n+ i| 2 ^sgn u n+ i (mod27r ) , 
written in the perception, 


7T 


£\ 3 / 2 37 r / a \ 3 / 2 71*7210 2/5 

' , Q = yr ( — ) Nh, e b = a \^N 2 b 


£b 


J 


4 \£ b j 


(58) 


(59) 


where u n+ 1 , 0 n+1 and u n , 0„, are respectively the values of the norined momentum and 
phase at times (n + 1)£ and nf. 

In writing (58) we have used the relationship between fields A w and B w , 

A/m = oh/ 2, b = B w /B 1 (60) 

which follows from equations (3). 

The quantity A/m has a clear physical meaning: it is the dimensionless representation 
of the ratio of the wave field work on one wavelength to the particle rest energy. 

Now we describe a relativistic particle motion in a time-like wave packet (TLP). 

In this case the Hamiltonian of the problem becomes 


H(p, z; 1, 9; t) = H 0 (p, I) + 


H a 


-A cos 0 5(t — n), 


(61) 


where we used the TLR of wave packet given by (25). Here 8(r—n ) is the Dirac function, 
r = [t/T], and T is the characteristic timescale. 

The associated with (61) equations of motion, written on framed manifold with the 
preceeding notation are 

p = k — 7 ^' B — Asin0 y ) 5(t — n), I = ■ A sin ^ ^ 5 (r — n), (62) 

Ha 


H 0 

z = pHf 1 , 9 = UBmHQ 1 . 
-0 = kz + 9 — uj. 


(63) 

(64) 


We retain in (63) only the terms of leading order. 

These equations can be transformed into a map by identification of the planes nT 
and ( n + 1 )T as follows. 

At first, one must note that the invariant of motion like that (53), lowering the 
number of degrees of freedom, is retained. Applying the invariant to (62 - 64)), we 
integrated these equations to find the coupled pair of nonlinear difference equations, 


p n+ i/m-p n /m+ (Nab/ 2) sin 0 n , N = [wT], 


07i+i = 077 + 


m 


\Pn+l\ 


sgnp n + i(mod 2-jt), 


where the variables p n , 0„ are taken at t = nT. 


(65) 
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In the representation 

f e \ 2 7 xa 2 bN 

U) =“■ Q = 


7 r 


2e? ’ 


where 

z b = (iV 2 &) 1/3 . 

equations (65) go over into a map, G n , 


( 66 ) 

(67) 


“1“ Q sin 'Ipn j 

G n : (68) 

f/Wl = Vti + 2(7r 3/2 /Q)|«n+ir 1/2 S5'n«n+l(mod 27 t), 
that is suitable for sequential analysis. 

The closed set of equations (68) is a measure-preserving map expressed in the terms of 
the canonical pair of variables ip, u. It is obvious, g n as well as G n inherit the canonical 
structure of the basic equations. Note that the maps g n ,G n are formalyzed by the 
essentially different transformations (59) and (66). In the following we refer to (58) and 
(68) as corresponding to the SLR case and the TLR case, respectively. 


5. Stochastic attractors of the system 

The maps g n and G n , or rather the families of maps depending upon the parameter 
Q, become suitable for sequential analysis. So thereafter we deal with high-frequency 
heating relativistic particles, expressed in terms of the discrete groups of automorphisms 
g n and G n , that act on smooth manifold. Denote by 

J = (69) 

the Jacobi matrix of these maps. It is important to note that the Jacobian of (58) is 
equal to one, therefore, g n and G n have a structure of the differentiable area-preserving 
maps. It stands to reason that u and ip are the canonical pair of variables. 

We now study the behavior of this dynamic system by computational analysis. We have 
numerically integrated equations (58) and (68) for several different values of Q from 
10“ 3 to 0.1. Figure 1 shows some our results computed for map (58) after (10 6 -10 4 ) 
iterations for one trajectory in the (u,ip) phase space. Shown in Figure 2 is the phase 
space of G n . The initial conditions were chosen in a random fashion and corresponded 
to the region of small values of (it, ip). 

We visual impression of the figures indicates that the phase curves are likely 
stochastic orbits. So, while smooth curves, piecewise smooth curves, and curves with 
finite length all have topological dimension one, the phase curves constructed by iteration 
function systems g n and G n have seemingly dimension close to two. Of course, to clarify 
this issue, we need to establish that dynamics of the system is realized on stochastic 
(strange) attractor. Let us consider a pair ( M,g n ) and/or a pair ( M,G n ), where M 
be a smooth manifold, and g n and G n are the differentiable area-preserving maps, and 
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either of these pairs determine the phase flow on M. From (58) and (68) follow that 
these pairs are invariable under the transformation 

ip — Y —ip, u ->■ -it, (70) 

and the inversion of a point with respect to a circle ?/>(mod27r) G S. On that ground, 
we define the following equivalence rules: 

(-7r,u) ~ (7r,u), (ip, -fr) - 00, ?r), (71) 

where ~ stands for the equivalence sign. 

The identification of these points is indicated by the arrows in Figure 1. 

The equivalence rules (71) allow us to represent phase space of the system as the quotient 
space having the quotient topology of torus, T 2 = S x S, ip(mod2ir ) G S,u G S [23]. 

Now we discuss the local topology considering a Jacobi matrix given by (69). Denote 
through Ai and A 2 the eigenvalues of the matrix J. Recalling that local structure is 
determined by the proximate neighbourhoods of fixed points, we have from (69) 

detJ = A r A 2 = l, (72) 

trj = Ax + A 2 = 2 + (tt 5 ^ 5 ) 1 / 3 , (73) 

for the SLR case; and 


det J = 1, (74) 

trJ = 2 + 7r 3/2 M" 3/2 , (75) 


for the TLR case. Here det J and tr J denote the determinant and the trace of this 
matrix, respectively. 

As pointed out in [11] the condition 

| tr J | = 3 (76) 


corresponds to a topological modification of a phase space, and its validity implies that 
the manifold have topology of a hyperbolic torus. Thereupon from (76) subject to (73), 
and/or (75) we find 


Ai 


3+ y/h 


3- y/E 


(77) 


Since Ai and A 2 are the dynamic characteristics such that A 2 > 0, Ai > 1, and the ratio 
Ai/A 2 is an irrational number, therefore the map (ip n ,u n ) — g n (ipo,Uo), (ipo^o) is an 
initial point of phase curve, forms a stochastic phase flow with the mean rate of a loss 
of infomation, K, 


K = InAi, 


(78) 


K is the Kolmogorov entropy. 

It is clear, the phase flow (M, G n ) possesses just the same properties. 
Present conditions (73), (75) and (76) ensure that the relation 


u b | = 7T 


(79) 
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take place, and it determines the upper bound of {«}. 

Like that the set {u, if} is a compact, whose structural stability is determined by the 
fractal dimension, df 

dj = l — In Ai/ln A 2 , dj = 2. (80) 

We call any compact a probabilistic fractal if its topological dimension is less than df 
and K > 0 [12]. In our case both these conditions are valid. 

As seen in Figures 1 and 2 for a given values of Q the stochastic region extends to values 
of u predicted by above equation (79). We know, that u depends on the parameters of 
this problem as given by (73) and (75). By that condition (79) determines in itself an 
equivalence class in the Q-parametric space. 

Considering (73,59) and (53), from (79) two equations result 
/I \ 2 / 5 

£b = a ( -6A t 2 J , (81) 

£ b = Qa K/w) 2 N 2 bj ~ ia (fV 2 &) 1/3 , (82) 

which determine the upper values of the energy spectra for the SLR- and TLR cases, 
respectively. Note the dependence e b on the driving field b is weak enough. This agrees 
with both the numerical solutions and the results of qualitative analysis. 

Finally, it should be noted that the phase flows, (M, g n ) and (M, G n ) are structural 
stable and typical because dj and K be an invariants, and sup{e} smoothly depends on 
the driving amplitude b, which is the controlling parameter of the system. Phase flows 
with such properties are said to be the stochastic (strange) attractors. Certainly the 
flows (M, g n ) and (M, G n ) be such attractors as was to be proved. 

At last, the achieved fractal measure df = 2 infers that points of phase curve 
evenly fill all obtainable phase space, or in other words, all states of our dynamic 
system are equivalent. This issue is supported by simulation in Figure 3 on which 
is represented a joint probability density, p('0,tt). The following algorithm was used for 
computing p. We partition all the phase space on the identical cells with the mesh size 
AipAu = (27r/40)x(27r/30). Like that p(4>,u) is proportional to the number of phase 
points in the element of phase space, AifA u. Figure 3 reveals that diffusion in the phase 
angle is very fast. Significant phase angle diffusion occurs on timescales order of tens of 
T, where T is the step of one iteration. On these timescales, the change in u is smaller 
(the numerical calculations indicate the characteristic time for establishing the uniform 
distribution in u is propotional to T/Q 2 at Q < 1), therefore, the variable u is a slow 
varying coordinate on the strange attractor. Statistical aspect of the problem will be 
studied in more detail in Sec. 6. 

6. Particle stochastic heating 

Particle dynamics in random electromagnetic fields is known to be described by the 
quasilinear theory (QLT) [25]. The QLT approach, in particular, has been employed to 
deal with diffusion of electron in the turbulent field of whistler waves packet [22]. 
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In coherent electromagnetic fields, on the other hand, the particle dynamics is not 
described by this theories. The nature of the diffusion in this case is the stochastic 
dynamics of particles, when the motion along stochastic trajectories gives rise to the 
so-called deterministic diffusion [1], General mathematical and physical aspects of this 
problem have been discussed, for instance, in the reviews [24]; [10]. 

The purpose of the present paper is to investigate electron dynamics in coherent 
packet of whistler mode waves. We have employed here the method in which the 
perturbation of particles orbits are included into a Fokker-Planck-Kolmogorov (FPK) 
equation [1]. 

In the present section, we shall exploit the canonical Hamiltonian structure that 
has been developed above. Thus, the drift kinetic equation follows at once since 
the canonical structure of equations of motion (58) and (68) has been established. 
As shown above, the maps g n and G n act on any strange attractor therefore the 
distribution function (probability density) f (u; t) on the attractor obeys the Fokker - 
Planck - Kolmogorov (FPK) equation 


df(u; t ) = 1_9_ 0/ 

dt 2 du du' 

which holds if Q -C 1 (Lichtenberg, 1983). 

Here D is the conventional diffusion coefficient in phase space, 


(83) 


D =< (u n+ 1 - u n f > T 1 , 


(84) 


in which (u n +i—u n ) is substituted from (58) or (68), and <• > denotes the phase average, 
T is the timescale of maps g n or G n . The function f(u,t ) belongs to the space of all 
differentiable functions supported in [ — 7T, 7r] , and the functional 

f f(u,t)du= 1 (85) 

J — If 

is the condition of normalization. 

First by means of g n ( or G n ) we calculate by formula (84) the diffusion coefficient 

D = Q 2 /2T. (86) 


Then, making use of result (79) proved above, along with (83) and (86), we evaluate 
the characteristic time for redistribution u over the spectrum 

t d ~ ul/D = 2 T{ir/Q) 2 . (87) 

A objective of this study is to determine the time-independent distribution function 
and the rate of heating. This requires a solution of the FPK equation together with 
a normalization (85). Thus a solution of the FPK equation along with the boundary 
condition /(— tt) = /(• tt) in the limiting case t > t, d may be given in the form of the 
uniform distribution 


f{u) = (2tt) \ uGl, 1 = (— 7T, 7T) 
for independent random variable u with mean zero. 


( 88 ) 
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Now, it becomes relevant to determine how the system evolves in time at t < t d . 
We newly exploit FPK equation with f(u) and its derivative df/du vanishing at the 
boundary. We introduce the moment < u 2 >= dim 2 f (u) , multiply equation (83) by 

u 2 , and integrate the resulting equation over u to obtain 


d < u 2 > 
dt 


= D. 


(89) 


Equations (83-89) describe two physically different situations: an electron motion in the 
SL wave packet, and, that respectively, in the TL packet. 

We know, that the variables u and e are related by (59) and (66) . This allows us to 
attach all possible states of e a probabilistic measure, namely, the probability density, 
f(e,t), which is associated with f(u,t) via the measure-preserving point transformation 


/M) = f{u,t)(du/d(e)). 


(90) 


Representing our results we start with the SLR case. In the case the variables u and e 
are associated by relation (59). 

Then from equations (88) and along with the norming 


[ f(e,t)de = l 
Jo 

we derive the steady-state distribution 

f(s)d€ = {e G R+ |0 < e < e b }, (91) 

for the random variable e with the mean value < e >— 0.6e;,, and the relative standard 
deviation 


y/< s 2 > - < £ > 2 


= 0.4. 


(92) 


V< £ 2 > 

Above equations (91) and (92) describe the state density and level of fluctuations in the 
energetic spectrum, e b is the upper value of this spectrum given by (82). 

Then we used relations (59) in equations (89) and (87) to find the rate of diffusion, 


d<e z > 


D, D 


a 3 N 2 b 2 


t < t d , 


dt 7 16 T 

and the characteristic time of diffusion in energy, 


(93) 


t d = 16 T (AT/64& 2 ) 275 . (94) 

The distribution function f(u,t ) in the TLR case obeys also the FPK equation as 
consequence of the fact that the equations of motion, expressed in coordinates (it, ip) 
reflect the underlying canonical Hamiltonian structure of equations (68). 

Now we use the dynamical equations (68), written in representation (66) and the 
local transformation (90), to derive equations, governing diffusion in energetic space. 
Thus the energy spectrum of particle follows from formula (88) 


f{e)ds = 


~2 ede , 

e b 


(95) 
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where e b is given by (82). We have respect to equations (83, 84) and (87) to find the 
equations, describing evolution of the energy spectrum 


Ft ne ' t) = \i m rs f{£j) ’ 


O(e) 


a 4 N 2 b 2 
32 Te 2 ’ 


(96) 

(97) 


d < s 4 > 
dt 


a 4 N 2 b 2 
16 T ’ 


t < t d , 


(98) 


t d = \T(N/bf\ (99) 

Notice expressions (93) and (98) describe a stochastic process differs essentially from 
the Brownian-like diffusion. Furthermore solutions of equations (93) and (98) predict 
qualitatively different behavior of the system. Thus in the SLR case the particle energy 
grows with time as i 1 / 3 ; in another case the energy growth is proportional to t 1 / 4 . 

Now we describe the effects associated with the stochastic heating of particle. First 
we discuss the pitch angle distribution over particle energies. Denote by x p the pitch 
angle and let x 1S its complementary angle. Then on account of A w /m <C 1 we have 

tanx = — ~ — . (100) 

v t s t 

Taking account of invariant of motion (53) we derive the following dependence 

tan x= {s z /oi) 112 . (101) 

Thereupon correspondence (91) along with the measure-preserving transformation 
f(x) dx = f (s') de yield the time-independent pitch angle distribution /(%), namely, 

/(x)^ = 3(a/£ fc ) 3/2 (l+tan 2 x)tan 2 xdx, X € (~X&, Xb), (102) 

tan Xb = (s b /a) 1/2 , 


and its dependence on the particle energy 

/(xOO) =3- 


el^a 1 / 2 


(103) 


Here s b is given by (81). A function given by (102) is the concave symmetric function 
of x, having the maximum at x = tt/ 2 (x p = 0), its derivatives tend to infinity as 
X -> ±7t/2. This function describes the so-called U-like distribution. 

In another case under condition (38) we write down expression (100) simply as 


X ^ s/a, 


(104) 


xe(-x&,x&), Xb = s b /a , 


(105) 
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f( Y ) = ^ Xft, 

0, elsewhere, 


(106) 


for the random variable x with mean zero, and its dependence on the particle energy 

/«£))= < 4 - ( 107 ) 

£ b 

Now the equation governing diffusion in pitch angle may be derived via the measure- 
preserving transformation f(x,t) = /(e, t) de/dx and the relation x = e/a. For /(y, t) 
satisfying these conditions on [— Xb,Xb], we write down the following FPK equation 


p(x,t) = 
B(x) = 


N 2 b 2 


32 Tx 2 ' 


From these equations we find the timescale for scattering, 


(108) 

(109) 


t. 


= 32 T 


xt 

N 2 b 2 


T ( N/bf ,z , 


and the scattering rate, 


d 4 N 2 b 2 
dt <X 16T ’ 


t < t s . 


( 110 ) 


( 111 ) 


According as (99) and (110), we note the timescales for heating and scattering are 
comparable. 

Next we show that the stochastic heating is accompanied by a spatial diffusion 
of electrons across a homogeneous magnetic field. Indeed, because of the one-one 
correspondence r = £t/u B and the invariant of motion (53), the distributions of electrons 
over r and e are related uniquely by the relation f(r,t) = f(£,t)de/dr. It means that 
f(r,t) obeys the FPK equation 




AM = 


a A N 2 b 2 


( 112 ) 

(113) 


327V 2 u;|’ 

which follows from equation (96) at once. 

The function f(r,t) is supported on [r 0 ,r&], where r 0 = l/u) B and r b ~ £b/(^B, and D t is 
the coefficient of diffusion of electrons across a magnetic field. 

One should be noted, that a classical Brownian-like diffusion process predicts that 
the mean square displacement grows linearly with time, and its rate is given by 
Dff ~ u ei < r 2 > representing the Bohm-like diffusion coefficient {v ei is the frequency of 
collisions). 
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Note that D'f has the strong dependence on B (as B~~ 2 ), while the coefficient of diffusion 
given by (113) is non-sensitive to change in B, and it hinges on b and e as b 2 , e~ 2 . 
Equations (112 and (113) are correct, if the inequality r b /L <C 1 holds. In view of 
equation (82), this requiring can be written in the form r b /L ~ (w/w b )(£ 6 /1V), (<C 1). 
For typical parameters this condition is satisfied trivially. 

We have done the solution of this problem. By way of illustration, we apply the results 
of our inverstigation to the Jupiter magnetosphere. 

7. Application 

It is known that the Jovian magnetosphere contains high-energy electrons spiralling in 
strong Jupiter^ magnetic field. Thus the existence of Jupiter’s radiation belts were 
confirmed via synchrotron emission measurements by Allen et al [13] and Fischer et 
al [14]. Of late Bolton et al [15] have bedf reported about the presence of electrons at 
energies up to 50 MeV in Jupiter’s inner radiation belts. They have been suggested that 
high energy electrons require substantial acceleration by processes other than adiabatic 
radial diffusion. The mechanism may be local acceleration of electrons by whistler mode 
waves, which observed in the Jovian magnetosphere [16]. 

To check the performance of the method, at first the parameters of the problem 
are chosen to reflect a typical whistler wave propagating along the L p = 6 shell, where 
the electron cyclotron frequency u B — 3.3 • 10 5 s -1 , the electron plasma frequency 
u) p = 2.2 • 10 6 s -1 , and the wave frequency to — 1.5 ■ 10 5 s -1 . Thus the dispersion 
relation (8) corresponds, undoubtedly to the problem, and by (8) we find the phase 
velocity, v p h — 0.08 and the wavenumber k = 6 • 10~ 5 cm -1 The theory is based on 
three parameters, namely, a, N, b. First we evaluate a = 2v ph co B /aJ, a = 0.32. Typical 
magnitudes of whistler mode waves are in range of 10 7 G, consequently the parameter b 
is of order 10 -5 [17]. To evaluate the parameter N = kL for the space-like wave packet, 
we need to estimate the characteristic size of the domain of resonant interaction (the 
zone of intense wave activity). Thus the characteristic size of the region is of the same 
order of magnitude as the spacescale of inhomogeneity of a background magnetic field 
near the Io torus, we have accepted L = 5 ■ 10 8 cm. By that the value of N is about 
N — 3 • 10 4 . In the following we will use these values to obtain an estimate of some 
quantities. First by means of (82) we calculate the upper value of energetic spectrum, 
E b ~ 3.5 MeV. Considering that the transit time of the particle through the intense wave 
activity domain is T — L/v z , i.e. T ~ 1.67 • 1CT 2 s, we evaluate by (94) the timescale 
of stochastic diffusion over the energetic spectrum, t d ~ 8 — 9 hours. Now the effect of 
pitch angle scattering should be evaluated. From (102) at e b = 7 and a = 0.32 follows 
that in this case whistler waves can diffuse electrons in a cone with the vertex angle, 
Xv = 2 Xb, which is about 156°. The equation (103) indicates a degree of scattering 
anisotropy, deg A = /(x)(e = £ b )/f(x)(£ = 1), increases with e up to deg A ~ 50. 

Now we apply our results to Jupiter’s inner radiation belts. Here we use the 
following values: co B = 2 • 10 7 s _1 , uj p = 5.65 • 10 5 s -1 , w = 1.5 • 10 5 s -1 . Thereafter we 
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obtain v P h — 0.95, and a, ~ 270. Now the transit time is L / v gr . where L is the lengthscale 
of inhomogeneity at these locations, L ~ 10 9 cm, therefore T is about 3- 10 -2 s. In view 
of (3), B w — E w /v ph , and using a spectral density of (1 — 3) • 10 _u V 2 m -2 /Hz obtained 
by Kurth et al [16] and Hobara et al [18] we find B w ~ 1.6 • 10~~ 8 G. Like that the 
reference parameters are: b = 1.6 • 10 -8 and N — [a >T] = 5 • 10 3 . Then we calculate by 
(98) and by (99) the upper value of the particle energy, Eb ~ 50 MeV, and respectively 
the timescale of establishing the energetic spectrum, td — 7 — 8 days. These expressions 
yield the mean rate of heating, namely, £&/2f rf = otb/T, Bb/2td ~ 100 eVs -1 . It is clear 
the mean heating rate depends linearly on the wave field magnitude, moreover the mean 
energy acquired by a particle equals to the wave field work on one wavelength multiplied 
by the ratio of the diffusion time to the characteristic particle transit time through the 
wave packet. We find that the phase space diffusion leads to a diffusion of electrons 
across an external magnetic field. Then (113) we calculate the rate the diffusion caused 
by chaotic motion, D t ~ 2 • 10 7 e -2 cm 2 /s. This process is more effective than a classical 
Bohm-like diffusion. To show this we use the following simple estimate. Taking plasma 
density of 10 2 cm -3 and assuming v ~ c, we obtain from D f ~ t/ e , < r 2 > that 
Df 3 ~ 5 • 10 -6 cm 2 /s. Thus the diffusion in phase space results in an enchancement of 
the absolute radial diffusion. Now we respect to the pitch angle scattering. The results 
(105) and (107) show that the wave packet to effectively scatter electrons in pitch angle 
leading to establishing the distribution peaked at Xp — 74° with the degree of anisotropy 
deg A = 100. Of particular importance is that in this energy range, the characteristic 
scattering time is comparable to the time of diffusion in energy. 

Relativistic charged particles lose their energy through the emission of electromagnetic 
waves. When the dominant loss channels are the synchrotron radiation and the inverse 
Compton effect, the rate with which a relativistic particle loses its energy is 

£ = -/3b 2 , (114) 

where j3 is the radiation coefficient. 

As a result, we arrive from (114) at the formula 

t r = {fie)- 1 = 5.1 • 10 8 (B 0 /B) 2 £-\ B 0 = 1.0 G, (115) 

for the characteristic lifetime of the particle [19]. 

Comparing t r and the timescale of diffusion over spectrum, G, shows that typically 
t r A> td, therefore the effect of emission on stochastic heating is negligibly small. Thus 
expressions (81) and (82) for upper values of particle energy spectra remain, and changes 
in the shape of spectra are very little. It is known [19] that wave spectrum of synchrotron 
emission is characterized by a well defined frequency, 

u r = wge 2 , ( 116 ) 

and it can be represented with the help of the relation (95) in the form of the uniform 
distribution 


f(u) r ) du) T = w b 1 du r , 


(117) 
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where co b = lub^I is the limiting frequency in the spectrum. 

We observe from (116) that the frequencies at which an electron emits the synchrotron 
radiation is considerable higher than the plasma frequency uj p . It is interesting in that 
this radiation, providing an information about the energy spectra, can leave the local 
plasma, and then be detected by an outside observer. 

Now to see that our approach correctly interprets the facts, it is necessary to find 
the range of validity of our results. The equations of motion (56) and (57) in the region of 
relativistic energies are obtained in the faithful representation, and the ultra-relativistic 
approximation, e 2 ;» max{a 2 , 1}, has been employed in deriving (58). According to 
(81), this inequality may be written in the form 

2 (uj B /u)v ph £~ l < 1. (118) 

Owing to (82) the latter can be transformed into the requiring 

b = B w /B > AN~ 2 oT 5 I 2 (~8-l(T 8 ). (119) 

In another case electron motion obeys equations (68), obtained in the approximation, 
1<£ 2 <q 2 . Taking into account expression (82), we write down this condition as 

8N~ 2 a~ 3 < b < 8N~ 2 . (120) 

By typical condition for Jupiter’s inner magnetosphere the above expression gives 
10“ 14 < b < 4 ■ 10~ 7 . These conditions are fairly easy to check. Recalling a spectral 
density measurements [16] and [18] yield B w ~ 10 -7 G, we conclude that (119) and 
(120) trivially satisfied, therefore, the upper limits for the energetic spectra given by 
(81) and (82) are correct. 

At last, the basic equations are applicable provided kL 1, and u) B T 1, and 
A w /m < C 1, which are typically fulfiled. 

It should be noted that the geometry of an ambient magnetic field does not play a role 
in the given problem, thus the lengthscale of an interaction region is typically smaller 
than the travelling path, i.e. an ambient magnetic field is always locally an uniform 
field. To this may be added the analytic model describing the motion of a particle out 
of the interaction region must include the radial diffusion of electrons conserving the 
first and second adiabatic invariants. 

We conclude that for typical values of the wave field in the Io torus, significant 
diffusion occurs on timescales of the order of a few hours for electrons with energies up 
to 3.5 MeV. In the inner radiation belts stochastic acceleration of electrons by whistler 
waves could provide a viable mechanism for energezation of electrons from energies near 
1 MeV to 50 MeV over a period of a few days. The results obtained indicated the 
dependence of maximal electron energy e b on the magnitude of wave field b is rather 
weak, although e b increases with increasing b. 

The obtained energetic spectra have used for evaluating pitch angle distribution 
functions over different energies. Stochasticity in the pitch angle distribution of 
relativistic electrons spiralling in a strong magnetic field near Jupiter is appeared as 
sharp peak near 90°; in relatively weak magnetic field at L p — 6 electrons are mainly 
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scattered by whistler waves almost along the direction of an external magnetic field. 
The time to scatter a particle into the scattering cone is comparable to the time of 
diffusion energy. 

The electrons involved into a process of stochastic acceleration produce radiation via 
the synchrotron mechanism. Thus energetic electrons emitting observed hectometric 
[20] and centimetric [15] synchrotron radiation respectively in Jupiter’s outer- and inner 
radiation belts according to (116) and (117) must have energy comparable to 3.5 MeV 
and 50 MeV. Thereto the observation of synchrotron emissions is obviously evidence for 
in situ particle acceleration. This is a further confirmation that the model is correctly 
predicting the main features of the observations. 

Chaotic motion driven by wave packet is also responsible for enchanced particle diffusion 
across an external magnetic field. This result in reasonable agreement with the 
experimental data. 

8. Summary 

A canonical Hamiltonian approach has been employed to deal with the whistler wave- 
electron interaction and the stochastic heating of high-energy electrons in magnetized 
plasmas. The time-like and space-like wave packet representations have used in deriving 
the equations of motion for relativistic electrons. 

The irreversible dynamics on a strange attractor was studied via the FPK equations 
with a coefficient of diffusion calculated from the equations of motion. The FPK 
equations describe such macroscopic effects as the stochastic heating, pitch angle 
scattering and radial drift of plasma particles. The effect of stochasticity on spectrum 
of synchrotron emission was also evaluated. 

The mathamatical structure is realized in the Jovian radiation belt events as follows. 
The main feature of Jupiter’s radiation belts is its striking stability in a dynamic state 
far from equilibrium. 

On the other hand the results indicate phase flow of our dynamical system is structural 
stable strange attractor. Chaotic motion on the attractor gives rise to an irreversible 
process (the so-called deterministic diffusion), which leads actually to establishing an 
steady-state energy spectra and results in such important and easily observable effects as 
the stochastic heating and pitch angle scattering of plasma particles. Under conditions 
typical of this mechanism, the heating region is determined by the boundaries of the 
attractor and the heating rate is governed by the nature of the kinetics, which in 
turn depends on the canonical variables on the attractor. The same conditions impose 
limitations on the timescales of macroscopic effects and feasible extent of heating, so 
that understanding of this conditions is of great practical interest. 
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Figure 1. The phase space of the map g n . (a) One single trajectory of length 10 6 for 
Q = 0.017T. The trajectory was started from the point uq = 10 -3 , ip o = 10 4 . The 
identification of points is indicated by the arrows, (b) Close to zero, approximation 
g n loses its validity. Therefore in the vicinity of u = 0 we must use the solution of the 
original equations (56) and (57) expressed in terms of the variables of u, ip given by 
(59). The figure shows the boundary of the chaotic region is well approximated by the 
condition (79). The parameters are the the same as in Figure la. 
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Figure 2. The phase space of the map G n . A single trajectory of length 10 6 for 
Q — 10 _3 7r. The trajectory was started from the point (10 -4 , 10~ 3 ). 
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Figure 3c 

Figure 3. Joint distribution p(ip,u) computed numerically via G n at several different 
number of steps. The parameters are the same as in Figure 2. The Figure shows 
p{ip,u) at times (a) t = 10 2 ; (b) 10 4 ; (c) 10 5 , respectively. The diagram at t = 10 2 
shows uniform distribution over ip is already established. Figure 3c indicates uniform 
steady-state distribution, implying that the dimension of the attractor is 2. 



